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^ ■ Abstract 



^y-j , Atom-molecule equilibrium for molecular formation processes is discussed for boson-fermion, 

§ \ fermion-fermion, and boson-boson mixtures of ultracold atomic gases in the framework of quasi- 

■ chemical equilibrium theory. After presentation of the general formulation, zero-temperature phase 
diagrams of the atom-molecule equilibrium states are calculated analytically; molecular, mixed, and 

^ ! dissociated phases are shown to appear for the change of the binding energy of the molecules. The 

^ ' temperature dependences of the atom or molecule densities are calculated numerically, and finite- 

■ temperature phase structures are obtained of the atom- molecule equilibrium in the mixtures. The 
0^ , transition temperatures of the atom or molecule Bose-Einstein condensations are also evaluated 



from these results. Quantum-statistical deviations of the law of mass action in atom-molecule equi- 
' librium, which should be satisfied in mixtures of classical Maxwell-Boltzmann gases, are calculated. 



. and the difference in the different types of quantum-statistical effects is clarified. Mean-field cal- 

culations with interparticle interactions (atom-atom, atom-molecule, and molecule-molecule) are 
^ ' formulated, where interaction effects are found to give the linear density-dependent term in the 



effective molecular binding energies. This method is applied to calculations of zero-temperature 
phase diagrams, where new phases with coexisting local-equilibrium states are shown to appear in 
the case of strongly repulsive interactions. 

PACS numbers: 03.75.Mn,05.30.-d,31.15.bt,82.60.Hc 



* mmut@tmu.ac.jp 

^ yabu@se.ritsumei.ac.jp 



1 



I. INTRODUCTION 



The experimental success of the Bose-Einstein condensation (BEC) {lj| of the trapped 
ultracold atomic gas has made much progress in the physics of quantum gases 0, 3, J], which 
includes Fermi-degenerate systems and Bose-Fermi mixtures. 

Recently, using the Feshbach-resonance method, molecular formations have been per 



and two bosons! 



formed experimentally in ultracold atomic gases for two fermionsjsl,!^,^] 
Bose-Einstein condensations have been observed for the thus created molecules of two 
fermions 



IG 



n 



11| and two bosons[9]. In these experiments, the molecular binding energies can 



be tuned by continuous changes of the applied magnetic fields; especially, bound molecular 
states can be changed into resonances by shifting the binding energies above the atom-atom 
scattering thresholds. 

One of the interesting applications of ultracold molecules is the observation of the BEC- 
BCS crossover in experiments with atomic Fermi gases: continuous crossover between the 
strong-coupling molecular BEC and weak-coupling BCS superconducting states 
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171]. In ultracold atomic-gas experiments, the crossover occurred by a change 



of the molecular binding energy through the Feshbach resonance method 
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l20(, and 



the experimental 



crossover 
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success has led to a lot of experimental and theoretical works on the 



231] and molecular BEC physics 2J, |25]. 



Molecules in optical lattice potentials are also an interesting problem. After the first 



observation of molecular formations using Rb[26l] , many experiments have been performed 
on the long-lived state of the ^^Rb molecule |27l] , ^'^K difermion 281] . and ^'^Rb-'^^K boson- 



fermion heteronuclear molecules 
^^Rb molecules 



Quantum degeneracy has also been observed in gases of 



301 ]. The phase structure of the lattice-trapped atomic gas, which is produced 



through interparticle correlations, has much interest in relation to the strongly correlated 
condensed-matter system described by the Hubbard model. Thus many experimental and 
theoretical works have been done: for example, on the superfluid-Mott insulator transition 



in bosonic 
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,l3^ 
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351 ] , fermionic 36|, l37|] , and boson-fermion systems 2 



Other theoretical studies on ultracold molecules include coherent molecular solitons 



39|, 



coherent pho toassociation 401]. and coherent intercondensate exchange between atoms and 



molecules 



In this paper, we discuss atom-molecule equilibrium for boson-fermion, fermion-fermion, 



2 



and boson-boson mixtures of ultracold atomic gases with quasichemical equilibrium theory, 



42 



43l ] for the quantum many- 
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3 ody 



an extension of classical chemical equilibrium theory 

problem, which was originally developed for the electron system in superconductors 
Our special interest is in applications of this method to boson-fermion and boson-boson 
mixtures, especially occurrences of atom or molecule BEG, as shown in some works for boson- 



fermion mixtures 
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471]. In addition, in contrast to the fermionic system, where 



many-body quantum calculations based on the microscopic model have been established, the 
singularities from the boson degrees of freedom sometimes cause problems in calculations of 
boson-fermion and boson-boson mixtures; quasichemical equilibrium theory can give definite 
solutions in such cases. As another interest of this approach, it should be pointed out 
that we can easily include interparticle interactions, especially atom-molecule and molecule- 
molecule ones, which are sometimes omitted in many-body calculations. The quasichemical 
approach gives the equilibrium structures in a less model-independent way. It turned out 
that the effects of these interactions should change the atom-molecule equilibrium structures 
drastically in the strong-coupling region. 

In Sec. II, a general formulation of quasichemical equilibrium theory is presented for 
molecular formation or dissociation processes in noninteracting atomic-gas mixtures, and in 
Sec. Ill, the method is applied to boson-fermion, fermion-fermion, and boson-boson mixtures 
and the atom-molecule equilibrium structures are shown at zero and finite temperatures. 
Special attention is paid to the condition on the occurrences of atom or molecule BEG; the 
shifts of the BEG transition temperatures are discussed from the molecular binding energy 
effects. In Sec. IV, the law of mass action, which is satisfied in chemical processes with classi- 
cal Boltzmann statistics, is shown to deviate in ultracold molecular formation or dissociation 
processes by quantum-statistical effects (the law of quantum mass action). In Sec. V, we 
extend the quasichemical theory to include interparticle interactions for s-wave scattering 
processes (three kinds of atom-atom ones, two kinds of atom-molecule ones, and one kind 
of molecule-molecule one in combination) in molecular formation or dissociation processes. 
In the mean-field approximation, the original six coupling constants of the interactions are 
shown to integrate into two parameters. It allows discussions of the interaction effects to be 
very clear. The formulations are applied to molecular formation or dissociation processes 
in interacting mixtures, and we discuss the change of the equilibrium structures at T = 
through interaction effects; new phases with coexisting local-equilibrium states are shown 
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to appear in the case of strongly repulsive interactions. Section VI is devoted to a summary 
and outlook. 

We should comment that a relativistic extension is also possible of quasichemical theory 
and applications to diquark condensates in quark matter have been done in j48 |. 



II. QUASICHEMICAL EQUILIBRIUM THEORY ON MOLECULAR FORMA- 
TION IN ATOMIC-GAS MIXTURES 



A. Molecular formation or dissociation process 

Let us take an atomic-gas mixture consisting of two atomic spices Ai and A2 with masses 
mi and 1712] in quantum statistics, 2 are bosons or fermions, so that we have three kinds 
of combinations: boson-boson (BB), fermion-fermion (FF), and boson- fermion (BF). 

To develop a quasichemical equilibrium theory, we consider a molecular formation or 
dissociation process in the mixture: 

Ai + A2< — > {A1A2) = M, (1) 

where M is a composite molecule with mass tum, which is bosonic in BB or FF mixtures 
and fermionic in BF mixtures. 

The mass defect of the molecule is defined as Am a/ = {tiim — fni — 7712). The bound- 
molecule {AniM < 0) is stable in both vacuum and gases, and has molecular binding energy 
AE = ArriMC^, where c is the velocity of light in vacuum. In contrast, the resonance 
(Amjv/ > 0) is unstable at least in vacuum; however, resonance states can exist stably in 
gases, so that we consider both bound-molecule and resonance cases. 

Here we take the quasiparticle picture wherein the system consists of atoms Ai and A2 



and molecule M, which are quasiparticles 491]. In this picture, two-body interactions between 
"bare" atoms bring about two-body correlations in the mixture; their major effects are the 
creation of the composite molecule M with binding energy AEm and atoms and molecules, 
which are quasiparticle in the mixture, and interact through residual interactions, which are 
generally different from the original interaction between bare atoms. These quasiparticle 
interactions are generally regarded to be weak, and we neglect them in the first part of this 
paper. In Sec. IV, we introduce the quasiparticle interactions and discuss their effects on 
atom-molecule equilibrium within the mean-field approximation. 



The equilibrium condition for the process ([T]) is given by 

yUl + /i2 - /iM = AEm, (2) 

where /ii,2,M are chemical potentials of Ai, A2, and M. 

The molecular binding energy AEm in (12]) is generally very small (~ io^(^'"^'') eV) in 
molecular formation or dissociation processes in ultracold atomic-gases, and it takes the 
same order of magnitude with the chemical potentials of atoms or molecules fia at ultralow 
temperature (T = yuK — nK). Thus, the term AEm cannot be omitted in ([2]). 

For free uniform gases, the particle densities are given by the Bose and Fermi statistics: 
1 r d^k 1 

1 f (i^k 1 

= (2^ y eie.-.^y^'BT + 1 = p^i^3/2(-/iJA;BT)(for fermion (4) 

where ks is the Boltzmann constant and Aj-.a is the thermal de Broglie wave length of 
particle Aa at temperature T: 

\ rriakBT 

The one-particle energy Sa in ([3]) and (j4]) is give by Ea = /i^k^/(2mQ,), where k is the 
wave-number vector of particle Aa- 

The Ba and Fa in and (jlj) are the Bose and Fermi functions: 

^ , , 1 x^^Mx . , 

^^^'^^ = T(A)Jo 
^ , , 1 /""^ x-^'^dx 

= rxiyio 

where z/ corresponds to the fugacity and T{A) is the gamma function. The Bose function 
can be written with the Appel function (f){z, s) as Ba{i^) = (piA, e'^), and the Fermi function 
is expressed with the Bose functions Fa{i^) = Ba{i^) — 2^^^BA{2h'). 

The Bose function Ba in ([6]) converges in z/ > (or /Xq, < in the chemical potential) 
and becomes singular at yU^ = 0, with which a phase transition occurs to the Bose- Einstein 
condensate of the boson A^- In the BEC region, the chemical potential vanishes and the 
condensed density n^°^ in ([3]) takes a finite value (it vanishes in the normal region). 

In the process ([1]), the particle-number conservations for 2 give the constraints 

ni+UM = ni^t, n2 + nM = n2,t, (8) 



where ni^t and n2,t are the total number densities of atoms 2, which consist of isolated 
atoms and constituents in the composite molecule M. 

Solving Eqs. ([2]) and ([H]) with and (jl]), we can determine the densities ni^2,M in 
equilibrium at temperature T for the parameters m^, ria^t = 1;2), and AE. 



B. Scaled variables 



We now introduce scaled dimensionless variables; they simplify the form of the above 
equations and greatly reduce the number of parameters. 

The scaled masses rha are defined by rha = nia/mM, where rriM = ttt^i + ^2 + Auim- In 
the ultracold atomic-gas system, the mass defect {ArriM ~ 10^^ — 10^^'' eV) is highly smaller 
than the atom or molecule masses (~ GeV); thus, we use the approximation itlm ~ rrii + m2 
(the conservation law of mass in chemical processes) throughout this paper. 

We introduce Ut = ni^t + n2^t and Eg = fi^iritf'/'^ /mM as scaling quantities for the particle- 
number densities and the energies, respectively: for example, ha = rta/nt (a = 1, 2, and 
M), T = ksT/Es, AEm = AEj^j/Eg, etc. Eg can be interpreted as having the meaning of 
the Fermi energy for the fermionic matter of fermions with mass tum with density rit, but 
we use it simply as a scaling quantity with the dimension of energy. 

Using the scaled quantities, the equilibrium conditions ([2]) and ([H]) become 

^1 + fi2 - fiM = AEm, (9) 
hi + hM = hi^t, h2 + hM = h2,t- (10) 



The Bose and Fermi statistics formulas and (jl]) become 

B^/2{-fijf) + (for boson A„), (11) 



rrini 



rriai 



ria = -^3/2(-/ia/T) (for fermiou A^). (12) 

Now the problem of obtaining the atom-molecule equilibrium states is reduced to solving 
Eqs. ([9]) and (fTOj) with (fTT!) and (fT2l) for temperature T, molecular binding energy AEm, 
and mass and total density of Ai (mi, hi^t)- The scaled atomic masses and the scaled total 
densities take the values of < rha,ha,t < 1 and satisfy 

mi +7712 = 1, hi^t + n2,t = l- (13) 



These constraints play the role of reducing the number of independent parameters. 



III. MOLECULAR FORMATIONS IN NONINTERACTING ATOMIC-GAS MIX- 
TURES 

A. Bose-Einstein condensation and Fermi degeneracy 

Before presenting the calculational results, we discuss two interesting quantum effects: 
BEC for bosons and Fermi degeneracy (FD) for fermions. 

As discussed in the previous section, the phase transition to the BEC of the boson 
occurs when /ia = in (fTTI) : the transition temperature Tc is 

~ 27r(n«)2/3 (n^)2/3 

Tc= ^ ' — C 3/2 -2/3 ~ 3.313^^^ — , (14 
nia rria 

where C(3/2) is the Riemann zeta function. The condensed and thermal parts of the number 
density at T < Tc are given by ( |TT|) : with scaled quantities, they become 

' 27r 



-f^=R^) C(3/2), 



1 - 



(15) 



From these equations, we can find that, if bosons Aa exist in the mixture at T = (n^ 7^ 0), 
all Aa should condense into the BEC (?2^°^ = at T = 0). 

Different from bosons, fermions go into the FD state at very low temperatures; especially, 
at T = 0, Eq. ( fT2l) becomes 

na = {l^F,a) ' , (16) 

where fip^a = f^aiT = 0) is the Fermi energy. 

The transition into the FD state is not a phase transition, but a continuous one, so 
that, different from the BEC Tc in f|T4l) . no clear boundaries exist for FD. Instead, as an 
estimation of the occurrence of FD, we use the temperature at jlF,a = 0: 

(67r2n^,)2/3 (n^)2/3 

Because of the permitted ranges of the scaled densities < n < 1, we can find Tc ~ Tp 
from Eqs. f|T^ and f|T71) as a rough estimation; quantum effects begin to appear at the 
same order of temperature in both noninteracting Bose and Fermi gases. 



B. Molecular formations in the BF mixture 

Let us consider the BF mixture {Ai = B, A2 = F) with the molecular formation process 
B + F ^ M = {BF). The atom or molecule densities nB,F,M in atom-molecule equilibrium 
are obtained from the equilibrium condition (jH]) with Eqs. fllOI) - fll2p . Interparticle interaction 
effects are neglected, which will be discussed in a later section. 

To understand the overall structure of the atom-molecule equilibrium, we consider the 
phase diagram at T = 0, which is obtained with analytical calculations as shown in Appendix 
B. Fig. [1] shows the T = phase diagram of the BF mixture with the same atom masses 
ffiB = rhp = 1/2 for the molecular binding energy {AEm) and a total density of B {nB,t)- 
Because of the constraints f[T^ . AEm and nB,t are unique parameters to determine the 
equilibrium state. The bracketed letters in the regions or lines in Fig. [1] represent the 
species which exist in equilibrium at T = 0; for example, {B, F, M) in the central region 
means that atoms B and F and molecule M coexist there, etc. 

1.0 
0.8 
„ 0.6 
0.4 
0.2 
0.0 

-20 20 

AEm 

FIG. 1: The T = phase diagram of BF mixtures with the same atom masses in the AEM-nB,t 
plane. The bracketed letters show what kinds of particles exist in equilibrium at T = 0. 

From Fig. [T|, we find that the BF mixture at T = shows the structure that, in the area 
left of the border BAC {AEm ~ 0), the states with molecules as many as possible become 
stable (molecular states) and, in the area right of the boundary BC, all molecules dissociate 
into atoms B and F (dissociated states). Between these two areas, mixed states of atoms 
and molecules become stable in the sense of equilibrium. 

In Fig. [H we can read off the existence condition of the BEC of atom B; it is in the 
regions with the bracketed B because the bosons always condense into a BEC at T = 0. 

We now turn to the atom-molecule equilibrium of the BF mixture at finite temperatures, 
which are obtained in numerical calculations using (jlj) and (Q. Fig. [2]shows the temperature 
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dependences of the scaled densities tib and hm in the BF mixture with masses = rhp = 
1/2 and the same total atomic-number densities nB,t = "^f,* = 1/2- The lines d, e, and / in 
Fig. [2] are for several values of the molecular binding energies: AEm = 0, —4.78, and 9.57, 
respectively. The BEC line f lT^ and the ftp = line, Eq. f[T7|) . are also drawn in Fig. [2] as 
lines a and b, respectively; they show the boundaries where quantum-statistical effects begin 
to occur. 

0.5 
0.4 
0.3 

m 

0.2 
0.1 
0.0 

-10 1 2-10 1 2 

logiof logiof 

FIG. 2: (Color online) Temperature dependences of the number densities ub of the bosonic atom 
B [(a), left] and um of the molecule M [(b), right] in the BF mixture with the same atomic masses 
and the same total atomic densities. The hues d (dashed line), e (double-dot-dashed line) and / 
(triple-dot-dashed line) are for mixtures with AEm = 0, 9.57 (resonance state), and —4.78 (bound 
state), respectively. The BEC Tq line of B (line o) and the fip = lines for F (line b) and M (line 
c) are also drawn. 

Fig. [2]shows that, at the high-temperature limit, n^/ converge to (complete dissociation 
of molecules into atoms). It can be explained by the energy-entropy balance in the free 
energy. The atom-molecule mixed states have a reduction in the energy part of the free 
energy; in contrast, the dissociated states have a reduction of the entropy part because the 
state density of -B -|- F is larger than that of molecules. With increasing temperature, the 
entropy contribution becomes large and the reaction process swings to the dissociation of 
molecules. This mechanism can also be applied to BB or FF mixtures because the binding 
energy and quantum-statistical properties are less effective at high temperature. 

The equilibrium states have a large dependence on the binding energy AEm around and 
below the BEC Tq and p^p = lines. At T = 0, they converge to states belonging to the 
phases on the horizontal 'hF,t = 1/2 line in Fig. [TJ the molecular, mixed and dissociated 
phases divided by the points, AEm = Ai = -{3n^f/''^/2 ^ -4.78 (point A) and A2 = 
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(Svr^)^/^ ~ 9.57 (point D). The temperature dependence in Fig. [2] is also classified with the 
values of AEm- (i) Ai?*/ < Ai (lines /). With decreasing T, the ub becomes small and 
goes below the BEC Tq line (no BEC of atom B for all temperatures), (ii) Ai < AEm < A2 
(line d). With decreasing T, ub (and also np) and hm run into the BEC and fl > regions, 
and converge into finite values, (iii) A2 < AEm (line e). With decreasing T, um takes the 
maximum value on the right of the fl = line (line c) and decreases to nM = along the 
line from below. 

The atom-molecule equilibrium in the BP mixtures can be explained by competition 
between the quantum-statistical effects and the binding energy of M; the molecule states 
give a free-energy reduction because of AEm < 0, but at low T, the molecules constitute 
the FD states, which have large kinetic energies coming from the occupied high-energy one- 
particle states. In contrast, in dissociated states, in spite of the kinetic-energy contribution 
from the FD states of the fermions E, the bosons B can reduce the kinetic energy largely 
as they condense into the BEC at low-T. Thus, depending on the positive or negative value 
of AEm, the dissociated or molecular state becomes stable, and a mixed phase appears 
between them. 

C. Molecular formations in the FF mixture 

In the FF mixture {Ai = El, A2 = E2), we consider atom-molecule equilibrium through 
the process Fl +E2 ^ M = {E1E2) in the same way as in the BF mixture. In this section, 
we take the FF mixture with the same atomic masses rfiFi = mpi = 1/2. 

Fig. [3] shows the T = phase diagram of the FF mixture. The notation for each phase 
is the same as that of the BF mixture (see Fig. [1]). Analytical expressions of the phase- 
boundary lines and points in Fig. [3] are given in Appendix B. 

The topological structure of the equilibrium phases in Fig. [3] is the same as that in the 
BF mixture (Fig. [1]); it consists of a "molecular phase" to the left of the border BAC, a 
"dissociated phase" to the right of the boundary BC, and a "mixed phase" between them. 

The existence condition of BEC of the bosonic molecules M is also read off in Fig. [3] as 
the regions with the symbol M in brackets. 

The equilibrium states at T 7^ can be expressed in the same way as in the BF mixture. 
In Fig. m we show the temperature dependences of npi and tim in the FF mixture with 
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FIG. 3: The T = phase diagram of the FF mixture with the same atom masses in the AEM-npi^t 
plane. The bracketed letters show what kinds of particles exist in equilibrium at T = 0. 




FIG. 4: (Color online) Temperature dependences of the number densities hpi of the atom Fl [(a), 
left] and njv/ of the molecule M [(b), right] in the FF mixture with the same atomic masses and 
the same total atomic densities. The lines c (dashed line) and d (double-dot-dashed line) are for 
mixtures with AEm = and 19.1 (resonance state), respectively. The jlp = lines and the BEG 
Tc lines are drawn for Fl and F2 (line a) and M (line b). 



f>'Fi,t = ^F2,t = 1/2 (the states on the line fipi t = 1/2 in Fig. [3D and AEm = and 19.1. 
In the FF mixture, the T = phases are classified with the value of AEm'- (i) ^Em < 
(molecule phase), (ii) < AEm < A3 (mixed phase), and (iii) A3 < AEm (dissociated 
phase), where the border points with AEm = and A3 = 2(37r^)^/^ ^ 19.1 correspond to 
point A and the crossing point between line BC and fipi t = 1/2 in Fig. [31 The two lines c 
and d in Fig. [H are for the border between two phases. 

We should comment on the relation _of^the present approach to the FF mixture with 
the BEC-BCS crossover theory 



21 



22 



23| . In the crossover theory, two kinds of bare 



fermions become dressed quasifermions, and quasimolecule states (or Cooper-pair states) 
appear as physical degrees of freedoms. The strength change of the attraction between bare 
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fermions gives the crossover between the BCS states (weak interaction) and the molecular 
BEC states (strong interaction). In the present quasiequilibrium approach, the fermions 
and the molecule should be considered as quasi-particles; the effects of bare-particle inter- 
actions are included as the existence of the molecule and its binding energy AEm- Really, 
the present approach is proved to give an approximated result on the strong interaction 
(BEC) side {AEm ~ 0) [2l|. In the weak interaction (BCS) side, the BEC Tc in the present 
approach vanishes at AEm = A^ as shown in Fig. [3] and disagrees with the BCS Tc- This is 
because the molecule states are treated as structureless bosons and no statistical correlations 
are introduced. It can be said that the present approach should give a good approximation 
in the strong interaction side and in the high-T region. 

D. Molecular formations in the BB mixture 

Before discussing the results, we should give some explanations of the singular property 
of the boson chemical potential at the double limit of ub = and T = in the no BEC 
region. 

Around the point of {ub.T) ~ (0,0), u = —^B/ksT becomes large, so that the Bose 
function S3/2 in the density formula ([6]) can be approximated by an asymptotic formula 
[Eq. flA9p in Appendix A]: then the chemical potential becomes 

^^B = kBT\Yl{XlnB) (18) 

with At oc T'^/^ defined in (El). Eq. ([I8]) shows that ^B is singular at {hb.T) ~ (0,0) and 
can take any negative value /iq if we take the limit through the pass ub ~ {\t,b)~'^&'^°^^'^^ ■ 
Now let us go ahead to the atom-molecule equilibrium in the BB mixture through the 
process Bl + B2 ^ M = {B1B2); then, the equilibrium condition ([9]) becomes 

fi-Bi + /iS2 - fj'M = AEm- (19) 
If all bosons condense into the BEC, the chemical potentials vanish, ^bi,B2,m = 0, at T = 0; 



it gives a contradiction in (fT9l) except for the case of AEm = (no triple BEC theorem) [48 1 . 
The negativeness of the chemical potentials determines the solutions of f|T9|) at T = 0: (i) 
AEm < (molecular phase): jiM = ^ (the BEC of M) and jla = AEm with = (a = Bl 
or B2). (ii) AEm > (dissociated phase): /isi = /iB2 = (the BECs of Bl and B2) and 
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fiM = —^Em with um = 0. The results are summarized in the T = phase diagram in 
Fig. [51 We should find that no triple BEC states exist except for the states on the boundary 
line of AEm = 0. Accordingly no mixed phase exists in the BB mixture. 
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FIG. 5: The T = phase diagram of BB mixtures in the AEM-nBi,t plane. The bracketed letters 
show what kinds of particles exist in equilibrium at T = 0. 
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FIG. 6: (Color online) Temperature dependences of the number densities n^i of the atom Bl [(a), 
left] and um of the molecule M [(b), right] in the BB mixture with the same atomic masses and 
the same total atomic densities. The lines c (short-dashed line), d (dotted line), e (dot-dashed line) 
and / (double-dot-dashed line) are for mixtures with AEm = +0, —0, and —1, respectively, 
where AEm = ±0 = hm,^+o ±e- The BEC Tq lines of Bl and B2 (line a) and M (hne b) are also 
drawn. 

Fig. [H] shows the temperature dependences of n^i and um in the BB mixture with rfiBi = 
^B2 = 1/2 and fisi^t = nB2,t = 1/2 (the states on the line riBi^t = 1/2 in Fig. E]) for 
AEm = ±0, ±1. The above-mentioned exclusive behaviors at T = for BEC have influences 
also at T 7^ as shown in Fig. [6j in AEm < 0, only the molecules become the BEC (lines 
e); otherwise, only atomic BECs can occur (lines /). The mixture with AEm = has a 
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singularity at point C; if we take the different kinds of limit AEm = ±0, BEC of atoms or 
molecules occurs, respectively. 

The exclusiveness between atom and molecule BECs in BB mixtures has been observed in 
a Cs experiment by the Innsbruck group j^; a sudden interchange occurs between the atom 
and molecule BECs when the resonance becomes the bound molecule using the Feshbach 
resonance method. This observation shows that the exclusiveness is incomplete and some 
mixture of atom and molecule BECs exists on the bound molecule side; it might be explained 
by an interparticle interaction effect (see Sec. V) or nonequilibrium effect. 



IV. LAW OF QUANTUM MASS ACTION 

In mixtures of classical Maxwell-Boltzmann gases, the molecular formation or dissociation 
processes through ([1]) satisfy the "law of mass action" : 

^ = K{T), (20) 

where ni,2,M are the number densities of the particle Ai^2,m and K(T) is an equilibrium 
constant, which depends on T, but not on rii t and n2^f 

In the present calculations of ultracold molecular formation or dissociation processes, 
quantum-statistical effects play an important role, so that they give the density dependence 
of K{T) in fl20l) : we call it the "law of quantum mass action". 

In order to make a comparison with the quantum cases, we derive an exact form of (120|) 
in the case of ideal Maxwell-Boltzmann gases, where the density of the particle Aa is given 
by 

= X-^ e^"/^^^, (21) 
with the thermal de Broglie wave length defined in ([5]) . Substituting Eg . (12T1) into the 



equilibrium condition ([2]), we obtain the classical law of mass action 



49|: 



(MB) (MB) . ,3 



T,M \ AG/keT 



K{T), (22) 



where AG is called the standard chemical affinity, and it becomes AG = AEm (the binding 
energy of the molecule Am)- 

To show the deviations from the law of quantum action by quantum effects, we evaluate 
the ratio of nin2/nM calculated in the previous section (including the quantum effects) to 

14 



n 



MB^MB I^MB 



Tin 



K{T) given by ([22D: 
ti = — 



n 



(MB) 
M 



nin2 



riM ni*^^)nf ^) riM K{T) ' 
The ratio R generally depends on T and fii^t because of the constraint ( |T3l) . 
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FIG. 7: (Color online) The ratio R defined by ([231) for the FF [(a), top left], BF [(b), top right], and 
BB [(c) and (d), bottom left and right] mixtures with the same atom masses. The mass defects are 
Attim = — 1 for the FF and BF mixtures and Aitim = —1 (c), +1 (d) for the BB mixtures. The 
lines are for temperatures: log^gT = oo (Maxwell-Boltzmann limit, solid line), 0.4 (long-dashed 
line), (short-dashed line), and —0.4 (dotted line), respectively. 

In Fig. [71 we show the ni j dependences of logi? for several temperatures (logT = 
— 0.4, 0, 0.4, oo) in the case of FF, BF, and BB mixtures with rhi = m2 = 1/2. The mass 
defects are Attim = —1 for the FF and BF mixtures and Attim = +1 [Fig. [71(c)] and —1 
[Fig. dZU^d)] for the BB mixtures. In the BB mixture, the behaviors of R are found to be 
very different for molecules and resonances. 

To understand the deviations from the law of mass action, we expand the ratio R in the 
high-T region analytically: 



R 



BB 



(24) 
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Rff ~ 1 + ^ (+At,i"'1 + At,2"'2 + >?t,m^m) , (25) 
Rbf ~ 1 + ^ (-Af. + A|,2'^2 - A| j^^n^/) , (26) 

where Rbb,ff,bf are for the BB,FF,BF mixtures. The leading-order terms i? ~ 1 in 
Eqs. ( p^ - (l26|) correspond to the Maxwell-Boltzmann hmit (law of mass action). Completely 
different contributions are found in the next-order terms (oc T~^/^) for Rbb and Rff in BB 
and FF mixtures, which are consistent with the results in Fig. [7](a) and Figs. [7](c) and El^d) 
at high T. The relatively small T dependence of the ratio Rbf for the BF mixture [Fig. [7](b)] 
at high T can be explained quantitatively from the cancellation of the boson and fermion 
contributions in the T^/^ term in (|26|) . 

V. EFFECTS OF INTERPARTICLE INTERACTIONS ON EQUILIBRIUM 
A. Interaction effects in the mean-field approximation 

In the previous sections, we discussed atom-molecule equilibrium in noninteracting atom- 
gas mixtures. If interactions exist between atom-atom, atom-molecule, and molecule- 
molecule (interacting) mixtures, they should modify the equilibration and sufficiently strong 
interactions could change the phase structures of the mixtures. 

Theoretically, the effects of the interactions for the free energy can be divided into two 
kinds: mean-field and correlation ones. The mean-field effect can be evaluated, for exam- 
ple, using the Hartree-Fock approximation, and it can be represented as interaction terms 
(background energy) in the one-particle energies of atoms or molecules in the mixtures. 
The correlation effects are defined as the contributions that cannot be introduced in the 
mean-field approximation. It should be noted that this division of interaction effects is a 
theoretical one and sometimes ambiguous. For example, in the BCS theory of supercon- 
ductors, if we take the normal quasielectron theory of a normal Fermi-degenerate vacuum 
in the Hartree-Fock sense, then Cooper pairs and their condensations are created by cor- 
relation effects; however, if we take the Bogoliubov quasi-particles as dynamical degrees of 
freedom, the BCS states can be understood as a kind of mean-field theory. In this paper, we 
discuss interaction effects for atom-molecule equilibrium in the Hartree-Fock-like mean-field 
approximation. 
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For interparticle interactions, we take the ones coming through the two-body s-wave 
scattering processes, which give dominant contributions in ultracold atomic gases except 
the strongly interacting or spin-degenerate ones. They can be introduced through effective 
interactions, for which we use the pseudopotentials 501] 



Vi,j = g^, E ^'(ra - r;^), {^.J = ^1, ^2, M) (27) 

between the ath and (3th. particles of the i and j species, respectively. The coupling constants 
Qij are determined from the s-wave scattering lengths between i and j species: 

9i,j = (28) 

where fiij = mimj/irrii + nij) is the reduced mass. 

Let us consider the contributions of the interactions in mixtures with equilibrium Al + 
A2 ^ M in the mean-field approximation. 

The interaction effects are introduced into the free energy as the background energy: 

F = Eo + E'"'- E (29) 

i=Al,A2,M 

where the Eq is the kinetic energy, which exists also in noninteracting cases. 

The contribution of the potential Vij in (1271) to the background energy i?™* is evaluated 
by i?*"* = (Vij). In the mean-field approximation, it is expressed by the number densities 
rii and rij. In the interaction between the same kinds of particles, it becomes 

EsB = ^-f {2n| - (ng))^} , Epp = 0, (30) 

for the boson B and the fermion F. The ng"* is the condensed density of B; when T ~ 0, 
we can use the approximation that n^jP ~ ub- The vanishing background energy for the 
fermion F originates in forbidden s-wave scatterings by Pauli blocking effects at T ~ 0. 
To avoid unnecessary complexity in the formulation, we redefine Qpp = for the fermions. 
The background energies coming from the interactions between different kinds of particles 
become 

Eij = QijUiUj {i ^ j). (31) 

From the above considerations, the total background energy in the equilibrium through 
Al + A2 ^ M can be approximated by 

^tnt ^ ^(my + gAl,A2nAinA2+ 9Al,MnAinM + 9A2,M^A2nM- (32) 
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,{0) 



Differentiating the free energy F in (129|) with respect to the density n^, we obtain the 
single-particle energy 

(33) 

j=Al,A2,M 

where ef^ is the kinetic energy of particle i: 



With the effective chemical potential defined by 



(34) 



Eq. (!33|) becomes 



j=Al,A2,M 



(0) / 



(35) 



(36) 



It should be noticed that the in f l36l) has the same form as that in the noninteracting case, 
so that the same Bose and Fermi statistic formulas (|TT1) and ( fT2l) can be applied also in the 
interacting case using the effective chemical potential /i^ instead of fi'f 

3/2 

n„ = — — B^i2{-ii'a/T) (for boson A^), 



2ti 
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^3/2 (-/u'a/^) (for fermion A^). 



Using eq. (l35ll . the equilibrium condition ([9]) for the interacting mixture becomes 



(37) 
(38) 



(39) 



where the effective binding energy AE']^,j is given by the density of the M molecule 

AE'j^ = auM + 7j 
where interaction effects are included in the two parameters a and 7: 



(40) 



tt= 51 gi,i + '^i9Al,A'i ~ 9ai,m ~ 9A2,m) 
i=Al,A2,M 



(41) 



7 = + {gAl,M - 9a1,AI - 9Al,A-2)'^Al,t + {.9A2,M - 9A2,A2 - 9Al,A2)'^A2,t- (42) 

In the derivation of eqs. (HU]1 - (H2|) . we have used the constraints (fTSjl . 
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As has been done in the case of noninteracting mixtures in Sec. IIB, we introduce scaled 
variables for the coupling constants and the effective binding energies: 

~9.^ = (43) 

^K: = (44) 

where rit = rii^t + n2,t and Eg = h'^{nt)'^^^ /tum- For the effective binding energy A£'^^, 
Eq. (HOl) becomes 

AE'^j = afiM + 7, (45) 
where the scaled parameters a and 7 are defined by 

«= E 9i,i + '^i9Al,A2- 9ai,m '-9a2,m)i (46) 
i=Al,A2,M 

7 = AE + {gAl,M - 9Al,Ai - 9Al,A2)^Al,t + {9A2,M - gA2,A2 - 9Al,A2)nA2,t- (47) 

Using the scaled variables, the atom-molecule equilibrium condition ( l39l) becomes 

fl2- f^M = ^E'm = OiflM + 7, (48) 

where Eq. (H5|) has been used. 

The atom-molecule equilibrium states of the interacting mixtures can be obtained from 
eqs. ( I37|) - (l39|) with the constraints (fT3|) . The interaction effects are included in the effective 
binding energy AE'j^^ in ( l40l) through the two parameters a and 7, which are determined 
from the coupling constants gi,j. Thus, we find that, in the mean-field approximation, one 
extra parameter is necessary in the equilibrium theory of interacting mixtures in comparison 
with that of noninteracting ones. 

B. Phase structure changes by interaction effects 

Before we give the numerical results for the atom-molecule equilibrium states of the 
interacting mixtures, we discuss qualitatively how the interactions shift and change the 
phase structures (PSs). We assume that all coupling constants are positive gij > to avoid 



possible instabilities from the spatial fluctuations of densities 38]. 

Different from noninteracting cases where the equilibrium condition (j9]) has a unique 
solution, the existence of the density- dependent term afiM in (j^8|) results in two or more 
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different solutions of (HHj) . which correspond to the different equihbrium states (coexisting 
phases). These solutions include locally stable and unstable states, so that we have to 
examine the behavior of Eq. (1481) and take out solutions corresponding to stable states. 

Close examination shows that two critical points > ad > i5c2 exist in the parameter a 
for the BF, FF, and BB mixtures, and the T = equilibrium structures can be classified 
into three regions with them as folloes, 

PSl (oci < Ci), Eq- (HSj) . has unique solutions for each value of the parameters a, 7, and 
77-1 2,t, and the T = phase diagrams become similar in structure with those for the 
noninteracting mixtures, including the dissociate, mixed, and molecule phases. (As 
shown later, the BB mixture is somewhat exceptional. ) 

PS2 (ac2 < a < ctci), Eq. fHHl) . has two stable (one unstable) solutions in the mixed phase. 

PS3 {a < 6lc2) has new coexisting phases appearing where both the mixed and molecular 
states become locally stable. 

Thus, interaction effects generally give complex phase structures in atom-molecule equi- 
librium with the phases with two locally stable states. The occurrence of these phases 
depends on the parameters a, 7, and 2,*, and the transitions of the states caused by the 
change of these parameters might be first order. In the phases with two locally stable states, 
the absolutely stable equilibrium state should be determined from a free-energy comparison 
of these states. However, the difference of the free energies of these states is generally small, 
so that effects that are not included in the present mean-field calculations (the correlation 
effect) may give comparative contributions, and, also, in real experiments, states which are 
not absolutely stable can occur through nonequilibrium and history effects. For these rea- 
sons, we should say that the stability of these phases is very subtle, and a study of the more 
detailed structures of these phases will not be done in this paper. 

From Eq. (H6l) . the negative contributions to a can be given by large values of ^ai,m 
and gA2,M', thus, atom- molecule interactions are found to be interesting and important in 
transitions to unstable phases. 
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FIG. 8: (Color online) The T = phase diagram of interacting BF mixtures with the same atom 
masses in the ^-ns^t plane for a = —20,0,20. The bracketed letters in the regions or lines show 
what kinds of particles exist in equilibrium at T = 0. 

C. Phase structures of the interacting BF mixture 

The critical values a^f and a^f '^^ for the BF mixture are given by 

1 + 



"c2 



2\ 2/3 
2\ 2/3 



1 



21/4^^/4 



1 + — 

TJlp 



(49) 
(50) 



Fig. [8] shows the T = phase diagram for interacting BF mixtures with = mp = 1/2 
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FIG. 9: (Color online) The T = phase diagram of interacting BF mixtures with the same atom 
masses in the ^-nB,t plane for a = —25. The bracketed letters in the regions or lines show what 
kinds of particles exist in equilibrium at T = 0. 
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for a = -20,0,20. Eqs. (09]) and ([50D give ajf^^ 20.7 and ajf^^ 28.7 when 

ifiB = rhp = 1/2, so that the phase diagrams in Fig. [8] are found to show the phase structure 
PSl; they have quahtatively similar structures as the diagram for the noninteracting mixture 
(Fig. [T]). Especially, the phase diagram for a = (solid line in Fig. [8]) is completely the same 
as that in Fig. [1] because of the vanishing density-dependent term afiM in ( HHI) . We also 
find that the boundary between the dissociate and mixed phases, CB, in Fig. [HI which is 
given by (]B17p in Appendix B, is independent of the values of a; this is because the um = 
condition on this boundary eliminates the a dependence in f l48|) . The boundaries between 
the mixed and molecular phases depend on a as shown in (1B15I) and (1B18I) in Appendix B. 
The a-dependence of the mixed phase can be understand from the position of the end point 
A in Fig. El 

which is given by flB20P in Appendix B. Eq. fl5T]) shows that the area of the mixed phase 
increases when a > and decreases in the case of a < 0. This behavior can be explained 
from fB5]) : in the case of a > 0, the molecule density hm has the effect of increasing AE'j^, 
which makes molecular formation difficult, and the area of the mixed phase becomes large 
(the a < case has the contrary effect). 

The phase diagrams with the phase structures PS2 and PS3 are given in Fig. [9] (a = —25) 
and Fig. [10] {a = —30). In Fig. [9l the occurrence of two locally stable states gives the 
additional two phases with the boundaries {AB + DF) and {DF + AB + EF) . The areas of 
these new phases in the PS2 structure are small; they are nothing more than substructures 
in the mixed phase. In the PS3 structure (Fig. [10]), the left-shifted phase boundary AB 
crosses the boundary BC, and new kinds of phase can be produced between AB and BC: 
region (b) in Fig. [TO] for example, where dissociated and mixed states coexist. Fig. [TO] also 
has a very complex substructure around the boundary BC. For extremely large values of 
a (Fig. [TT] a = —70), the regions with substructures shrink to the small areas around the 
vertexes B and C in Fig.[TT], and the whole phase structure becomes simple again; the phases 
where the dissociated and molecular states coexist appear in the central region, where the 
dissociated states include no molecule and the molecular states have as many molecules as 
possible. 
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D. Phase structures of interacting BB and FF mixtures 



In this subsection, we briefly sketch the structures of the T = phase diagrams of 
interacting BB and FF mixtures. 

The T = phase diagrams of interacting BB mixtures with m^i = 'mB2 = 1/2 are 
shown in Fig. [T2](a) {a = 20) and Fig. [T2r b) {a = —20). In the case of a > 0, the phase 
diagrams have the mixed phases [the region ABC in Fig. [TWa)]. which do not exist in the 
noninteracting cases. The position of the end point A is given by (]B25P in Appendix B: 
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FIG. 10: (Color online) The T = phase diagram of interacting BF mixtures with the same atom 
masses in the j-fiB^t plane for a = —30 [(a), top left], (b) and (c) (top right and bottom left) show 
the detailed structures of (a), (d) (bottom left) is the legend that explains the phases {i)-{n) in 
(a)-(c). 
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FIG. 11: (Color online) The T = phase diagram of interacting BF mixtures with the same atom 
masses in the 7-n_B,t plane for a = —70, where the phase (i) and (m) are coexisting phases of 
{B,F), {F,M) and {B,F), {B,M) as explained in the legend in Fig.[10];d). 
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FIG. 12: (Color online) The T = phase diagram of interacting BB mixtures with the same atom 
masses in the 7-n_Bi,t plane for a = 20 [(a), left] and a = —20 [(b), right], where the phases {i) and 
(m) are coexisting phases of {Bl,B2), {B2,M) and {Bl,B2), {Bl,M). 

the phase structure in a > is just PSl. When a < 0, Eq. (!52|) shows that point A moves 
across the boundary BC and locates to the right of it; the mixed phase disappears and new 
phases occur with coexisting locally stable equilibrium states [Fig. [TWb)]. That means that 
the phase structure is PS3 for a < 0. As a result, we find cicf^'' = c^if^^ = and no PS2 
phase structures exist in the interacting BB mixtures. 

In the case of interacting FF mixtures, the critical values a^i and become 

aff^^ = ("iFi^^ + (53) 
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air' = -2"^f^)'"^^, (54) 

which become ftcf^'* ~ —25.5 and q;|.^'^'* ~ —38.2 in the case of m^i = mp2 = 1/2- 

In Fig. [131 we show the T = phase diagrams for interacting FF mixtures with 
rhpi = = 1/2- The change of the phase structure is essentially similar to that in 
the interacting BF mixtures. In the PSl structure [Fig. [T^ a)]. the phase structures are 
obtained by deformation from the noninteracting ones. The substructures with the coex- 
isting states appear in the mixed phase in the PS2 structure [Fig. [T3r b)]. The end point 
A in Fig. [T3](a) and [T3l(b) crosses the boundary BC between the mixed and dissociated 
phases and new phases appear in the PS3 structure with the coexisting equilibrium states, 
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FIG. 13: (Color online) The T = phase diagram of interacting FF mixtures with the same atom 
masses in the ^-hpi^t plane for a = 20,0, —20 [(a), top left], —30 [(b), top right], —40 [(c), bottom 
left], and —70 [(d), bottom right], where the phases (i), (m), and (n) are coexisting phases of 
(F1,F2), (F2,M); (F1,F2), (F1,M); and (F1,F2), (Fl,F2,Af). 
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which have complex substructures [Fig. [T3](c)]: the structure becomes again simple for large 
negative values of a [Fig. [TST d)]. 

VI. SUMMARY AND OUTLOOK 

We have developed a quasichemical equilibrium theory for the molecular formation or 
dissociation processes in BF, FF, and BB mixtures of ultracold atomic gases and discussed 
atom-molecule equilibrium in these mixture. 

The law of mass actions has also been examined for the mixtures; it is satisfied well at 
high T. We have shown that the quantum-statistical effects of the atoms and molecules, 
which become more effective in ultracold temperature, give deviations from the law (law of 
quantum mass action). The quantum-statistical effects are shown to give different deviations 
for bosons and fermions, and, in BF mixtures, both contributions have a tendency to cancel 
out at high T. 

We have also discussed the effects of the interparticle interactions in the mixture within 
the mean-field approximation at T = and evaluated the shifts of the T = phase struc- 
tures of atom-molecule equilibrium in the mixtures. Especially, in the case of large repulsive 
interactions between atoms and molecules, the phase structures have been shown to change 
qualitatively with the occurrence of coexisting local-equilibrium states. We have given the 
conditions for the coupling constants with which the phase- structure changes occur. The 
atom-molecule equilibrium in interacting mixtures can be calculated also at finite tempera- 
tures within the present framework, and the results are planned to be published in another 
paper. 

The other correlation effects beyond the mean-field approximations should also be 
important-for example, in the BCS-BEC crossover problem. Combining the method of 



the Beth-Uhlenbeck approach [5l|, [52] with the present quasichemical equilibrium theory, 
we should discuss the correlation effects and the crossover problem from a less model- 
independent point of view. A study along these lines is now ongoing and will be presented 
in the near future. 

The authors thank T. Suzuki and T. Takayama for many useful discussions. 
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APPENDIX A: ASYMPTOTIC BEHAVIORS OF BOSE OR FERMI FUNC- 
TIONS 



In this appendix, we derive some formulas of the asymptotic behaviors of the Bose and 
Fermi functions (jH]) and ([7]) at z/ ~ and z/ ~ ±00. 

Before we discuss the asymptotic behavior, we prove the relation between Ba and Fa'- 

FA{iy) = BA{y) - 2'-^Ba{2u). (Al) 

This formula is obtained by integrating both sides of the equation 

= (-^2) 

Let us go to the asymptotic behavior of the functions Ba{i^) and Ba{i^) at 1/ ~ 00. Using 
the expansion of the integrand of Ba, 

„A-l 00 



X 

we obtain 



- = (A3) 

1 ; 1 



00 p — kv 00 p—kv 

-^-'^-"''-i:^. (A4) 



fOO 

k=i nA) Jo 

where we have used the formula of the gamma function: 



r{A) = / z^-^e~'dz. (A5) 
Jo 

Substituting flA4p into flAip . an expansion formula for -Fa(^) can be obtained: 

^a(^) = fli-l)'-''^. (A6) 
k=i 

Taking the first terms in Eqs. ( 1A4[) and ( ]A6p . we obtain the asymptotic behaviors of Ba{i^) 
and FaIi^) at ~ 00: 

Ba{v) ~ Fa{v) ~ e-r (A7) 

We turn to the asymptotic behavior of Ba{i^) around u = 0. When < a < 1, the 
point z/ = becomes an irregular singular point, so that only the asymptotic expansion is 
obtained for Ba- For this purpose, we consider the Hankel-type complex integral 

If z az ^^g^ 



r{A) J{oo;-u+) - 1 
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FIG. 14: Hankel-type integral path in the complex-number plane. 



The integration path (oo; — z/+) and the cut hne on the positive real axis are shown in 
Fig. [m The phase branches of z^~^ are fixed at x^~^ or e'^'^^''x^~^ on the upper or lower 
parts of the cut line. The integrand function has a pole at 2; = —u: 



6^+'" - 1 Z+U 



(A9) 



The integration path does not pass the singular point z = —u, so that we can evaluate 
the integral using the expansion ( 1A2I) : 

I = y ^— / z^-^e-^^dz = y ^(e^-^^ - 1), (AlO) 

,tl r(A) i(oc;0+) tx 

where we have used the Hankel-integral representation of the gamma function. 

On the other hand, we can deform the integration path and divide it into the path C-y 
circulating around the singular point z = —v and the Hankel-type path (00; 0+) around the 
origin O: 



z^~'^dz 



C-^ J{oo-fi+) 



(All) 



T{A) J{oo--u+) e-+^ - 1 T{A) 
The integral on C^i, is evaluated by the theorem of residue, and that on (00; 0+) can be 
attributed to the real integral: 

1 r z^~^dz 2TTi 



T{A) Jc-, e^+^ - 1 V{A) 



(A12) 



^ {e^^'^-l)BA{v). (A13) 



T{A) J(oo-fl+) e^+^ - 1 T{A) Jo 6^+*^ - 1 
Combining Eqs. flAlOl) . flA12p . and flA13l) . we obtain the asymptotic expansion of Ba(i^ 



Using the expansion of the series part in (lA14p . 

00 -ku 00 1 00 (_l.,,\n 00 (_,,\n 00 1 00 n 

EV = ^FiS^ = ^^Si^ = ^(-'r<(A - n)!^, (A15) 

k=l k=l n=0 n=0 "'■ k=l n=0 "'■ 
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we obtain the power expansion formula by Qpechowski[53l. |5J]: 



-^^ + E(-irc(A-.)^ 



(A16) 



Taking the leading term of the Opechowski formula, we obtain the asymptotic behavior of 
Ba{i') around u = 

' aA) (A>i), 

Ba{v) ~ < - Inz/ = 1), (A17) 

.ifefS + ^(^) (o<^<i)- 

Using Eq. ( lAll) . the asymptotic formula of FAiy) around i/ = becomes 

F^(z/)~(l-2i-^)C(A). 



(A18) 



It should be noted that the residue term in Ba{i^) is canceled out in ( lAlSP : this is consistent 
with the fact that no such singular terms exist in FaIi^) originally. 

The asyinptotic formula of -Fa(^) at z/ ~ — oo is obtained by the Sommerfeld expansion 



formula 
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(l)'{u)du 







+ 



-~0(a) + ^2F2„(O)0(2")| 

^ n=l 



a 



(A19) 



where is a cx)-differentiable function and the coefficients -F2n(0), the Fermi function at 
z/ = 0, are represented by the Bernoulli numbers B^- 

(1 -2i-2")(27r)2" 



2F2„(0) 



{2n)\ 



-B„ 



Using Eq. (lAlQ^ . the asymptotic formula of -Fa(z/) at z/ ~ — oo is obtained by 



r(A + i)' 



(A20) 



(A21) 



APPENDIX B: DERIVATIONS OF T = PHASE DIAGRAMS IN BE AND EE 
MIXTURES 

1. Noninteracting BE mixture 



Let us consider the mixed phase of the BP mixture of atom masses rfiB^F, which corre- 
sponds to the central region in the phase diagram (Fig. [I]). The equilibrium condition 
becomes 







2\ 2/3 



1 



V2 J niF 



2/3 



'3^2\2/3 

71 j 



= AE 



Ml 



(Bl) 
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where we have used fiB = (BEC state) and the Fermi energy formula (fT6|) for F and M 
with niM ~ 1. 

The boundaries AB, BC, and CA in Fig. [T] are obtained by setting tib = 0, tim = 0, and 
hp = 0, respectively: 

AB : AEm = ^i67r2(l - 2n^,,)]'/' " ^[6vrVt]'/', (B2) 
BC: A^M = ^[67rV,t]'/', (B3) 
CA: AEM = -^[67rV,i]'/'. (B4) 

and the end points of the boundaries, A, B, and C, are obtained from ( ]B4I) with n^^ j = 1/2, 
fIBSI) with fii?^ = 1, and fIBSI) with fi^^'t = 0: 

A : AEm = --^ 4.8, (B5) 

B : AEm = ^{Qrr' f' - ^, (B6) 
C : AEm = 0. (B7) 

The boundaries and end points in Fig. [1] are obtained from the above formulas for mp = 1/2. 



2. Noninteracting FF mixture 

In the case of the phase diagram of the FF mixture (Fig. [3]) with masses rfipi and mp2, 
the equilibrium condition in the mixed phase becomes 



37r 



2\2/3 , /o„2\2/3 . 

, J ^{hp^f'+r-^] ^fe)^/^ = A^M, (B8) 
\V2/ mpi \v2/ mp2 

where we have used fiM = ^ (BEC state) and the Fermi energy formula f|T6l) for Fl and F2.. 

The boundaries AB, BC, and CA in Fig. [3] are obtained by setting hp2 = 0, um = 0, 
and hpi = 0, respectively: 

[67r2(l - 2hp,,)f', (B9) 



AS : 




1 

2mi7i 


BC : 


AEm = 


1 

2rhpi 


CA : 


AEm = 


1 

2rhp2 



[QTT^hpi.tf' + -^[Qrr^hp2,t]^/\ (BIO) 
2mi?2 

[67r2(2n^i, - l)]2/3, (Bll) 
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and the points A, B and C are obtained from (1B9P with fipi^t = 1/2) (1B9P with ni?i^4 = 0, 
and fIBllI) with npi^t = 1: 



A : AEa/ = 0, (B12) 
B : AEm = ^{QrrY' - (B13) 

C : AEm = ^(67r2)2/3. ^ Zl^. (B14) 
2'mp2 mpi 

The boundaries and end points in Fig. [3] are obtained from the above formulas with mpi = 
mF2 = 1/2. 



3. Interacting BF mixture 

The T = phase diagram of the interacting BF mixture with masses rhs^F is obtained 
by the equihbrium condition (148|) for T = 0: 

In the case of a > a[i^\ the boundaries and end points of the phases can be obtained in 
the same manner as those in noninteracting cases. 
The results are 



2\ 2/3 



AB : 7 



'3f\ 
.V2 



rriF 



anB,t, (B16) 



BC: 7=(^)''':^4(*^ (B17) 



-2/3 

niF 

''?-^2\ 2/3 
OTT \ ^2/3 



CA: 7 = -(^_j (B18) 
for the phase boundaries {nB,t + ^F,t = 1), and 

^: 7 = -bl -2~-4.78-2, (B19) 



C : 7 = 0. (B21) 

The boundaries and end points in Fig. [7] are obtained from the above formulas with rfiF = 
rfiB = 1/2. 
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4. Interacting BB mixture 



The T = phase diagram of the interacting BB mixture with masses ^131,32 is obtained 
by the equihbrium condition fHS]) . 

At T = 0, all chemical potentials appearing in fl48l) . which are bosonic, can take two 
alternative possibilities: /t'^ = with Uk ^ (BEC) or /i^ < with = (/c = Bl, B2, M). 

In the mixed phase, the BEC conditions fi'^i = = I^'m = Si^^ afiM + 1 = 0; the 
boundaries, AB, BC, and CA, are obtained by substituting nB2 = 0, fiM = 0, and fiBi = 0, 
respectively: 

AB: 7 = -anBi,t, (B22) 
5C : 7 = 0, (B23) 
CA: j = -a{l-nBi,t), (B24) 

where we have used the constraints fiBi + nM = nBi,t and nB2 + = nB2,t- The end points 
are obtained by 

A: 7 = -f , (B25) 

B: 7 = 0, (B26) 

C : 7 = 0. (B27) 

It should be noted that the boundaries and the end points are independent of the atom 
masses rhBi,B2- 

5. Interacting FF mixture 

The T = phase diagram of the interacting FF mixture with masses m^.F is obtained 
by the equilibrium condition (l48ll . for T = 0: 

.2/3, ^(^^2/3 



'%Y'^inFir/'+(%Y'^inF2r' = anM + r (B28) 
V2 / rnpi \ V 2 / "^F2 



In the case of a > a).i , the boundaries and end points of the phases can be obtained in 
the same manner as those in noninteracting cases. 
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The results are 



AB : 


7 = 


1 

2m Fi 




(B29) 


BC : 


7 = 


1 

2171 PI 




(B30) 


CA: 


7 = 


1 

2mF2 


[67r2(2n^i,i-l)]2/3 + «(n^,,-l), 


(B31) 



for the phase boundaries (n^,* + = 1), and 

(B32) 

(6vr2)2/3 ^ Zli, (B33) 
(6vr2)2/3. ^ Zli. (B34) 

The boundaries and end points in Fig. [TT] are obtained from the above formulas with ifipi = 
mF2 = 1/2. 



A: 


7 = 


a 

~2' 


B : 


7 = 


1 

2miri 


C : 


7 = 


1 

2m F2 
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